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■ ABSTRACT 

\l ' HD 167971 is a triple system consisting of a 3.3-day eclipsing binary (05-8 V + 05-8 V) and an 08 supergiant. It is also a well 

known non- thermal radio emitter. We observed the radio emission of HD 167971 with the Very Large Array (VLA) and the Australia 
' Telescope Compact Array (ATCA). By combining these data with VLA archive observations we constructed a radio lightcurve cover- 

' ing a 20-yr time-range. We searched for, but failed to find, the 3.3-day spectroscopic period of the binary in the radio data. This could 

^ I be due to the absence of intrinsic synchrotron radiation at the coUiding-wind region between the two components of the eclipsing 

, binary, or due to the large amount of free-free absorption that blocks the synchrotron radiation. We are able to explain many of the 

observed characteristics of the radio data if the non-thermal emission is produced in a coUiding-wind region between the supergiant 
and the combined winds of the binary. Furthermore, if the system is gravitationally bound, the orbital motion occurs over a period of 
[ ~ 20 years, or longer, as suggested by the long-term variability in the radio data. We argue that the variability is due to the free-free 

, absorption that changes with orbital phase or may also in part be due to changes in separation, should the orbit be eccentric. 

^ . 

, Key words, stars: individual: HD 167971 - stars: early-type - stars: mass-loss - stars: winds, outflows - radio continuum: stars - 

■""^ ■ radiation mechanisms: non-thermal 



O 1- Introduction 

H . 

' Radio emission from most hot stars is due to thermal free-free 
emission by material in their stellar winds. A significant frac- 
tion of these stars, however, also show evidence of non-thermal 
. ^ emission. The radio fluxes of non-thermal sources are character- 
ized by a zero or negative spectral inde?([3 and a high brightness 
?H temperature. 

5^ For Wolf-Rayet stars, Dougherty & Williams (2mU) showed 
that non-thermal emission is strongly correlated with binarity. 
In a binary, the stellar winds of both components collide. On 
either side of the surface where the wind pressures balance, 
shocks are formed that define the extent of the colliding-wind 
region (Eichler & Usov I19931 I. Around these shocks, relativis- 
tic electrons are accelerated by the first-order Fermi mechanism 
(Bell |1978] l. These relativistic electrons then spiral in a magnetic 
field and thereby emit synchrotron radiation, which we detect as 
non-thermal radiation. 

For O stars, the link between non-thermal emission and bina- 
rity was, until recently, less clear. There are a number of spectro- 
scopically single O stars that are non-thermal emitters. In a sin- 
gle star, relativistic electrons could be accelerated in shocks that 
are due to the instability of the radiative driving mechanism (e.g. 
Owocki et al.; 19881 l. However, Van Loo et al. (I2006I I showed that 
this embedded shock model cannot explain the observed spectral 
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' The radio spectral index a describes the power-law behaviour of 
the flux: F^ocv"cc^-". 



index. They conclude that, just as for the Wolf-Rayet stars, all O- 
star non-thermal emitters should be binaries. Non-thermal emit- 
ters that are apparently single must therefore be binaries with 
an unfavourable inclination angle or a long-period orbit, mak- 
ing spectroscopic detection of their binary nature difficult. An 
example of such an object is HD 168112, which we studied in 
a previous paper (Blomme et al. 2005 hereafter Paper I). The 
periodic nature of the non-thermal radio fluxes suggests that this 
spectroscopically single star is in fact a binary, with a period es- 
timated to be 1 .4 yr. 

Problems remain, however, in understanding the observed 
radio emission in these colliding-wind binaries. Cyg OB2 No. 5, 
for example, is a 6.6-day period binary where a colliding-wind 
region is clearly seen (in high-resolution VLA observations) be- 
tween the binary and a third star The radio flux of the binary 
alternates between a "high" flux state and a "low" flux state and 
it does so on a time-scale of ~ 7 years, rather than the 6.6-day 
orbital period (Miralles et al. 119941 Contreras et al. '1997). To 
improve our understanding of the radio emission in binaries, a 
more detailed study of this type of object is therefore warranted. 

In this paper we study the non-thermal radio emitter 
HD 167971 (RA= 18*'18"05!895,Dec = -12°14'33'.'31, J2000). 
Leitherer et al. ([1987) found that this system consists of a 
3.3213-day period eclipsing binary and a third companion. The 
binary components are very similar to one another in mass and 
temperature, but the spectral types are not well determined: they 
can only be constrained to the range 05-08. From the abso- 
lute magnitudes, Leitherer et al. tentatively estimate both com- 
ponents to be main sequence. From a more detailed modelling 
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of the lightcurve, Davidge & Forbes (1 19881 1 propose a giant or 
supergiant classification. The spectral lines also show the third 
component to be a more luminous star, probably an 08 super- 
giant, and it is not known if it is gravitationally bound to the 
system or is just a line-of-sight object. No high spatial reso- 
lution observations are available, so the angular separation be- 
tween the binary and the third component is not known. Further 
orbital parameters of the binary are not well determined: the ec- 
centricity is probably e ^ 0. No detailed coverage of radial ve- 
locities as a function of phase in the binary orbit is available, but 
spectra taken near quadrature indicate a velocity amplitude of 
~300 km s"'. No significant velocity variations in the spectral 
lines of the third component have been detected. The distance 
to HD 167971 can be derived from its membership of the open 
cluster NGC 6604, which is part of the Ser OB2 association. 
From the distance modulus listed by Humphreys (1978) a dis- 
tance of 2 kpc is found. 

The purpose of this paper is to study the colliding-wind re- 
gions in the HD 167971 system, by combining new radio ob- 
servations with archive data. We construct a simple model that 
gives a plausible explanation for the important features seen 
in the radio fluxes and that is also compatible with the X- 
ray and VLBI (Very Long Baseline Interferometry) observa- 
tions. Detailed models, such as those developed by Dougherty 
et al. (I2003I I and Pittard et al. (I2006I I are not used here. Many of 
the orbital and stellar parameters are not known (especially for 
the interaction between the binary and the supergiant), making 
the parameter space too large for sophisticated modelling. 

The rest of this paper is structured as follows. In Sect. |2] we 
present the radio observations, with the details of the data reduc- 
tion split off" to Appendix lAl We discuss the variability detected 
in the radio in Sect. [3] The observations are interpreted in Sect.|4] 
and the conclusions are drawn in Sect. |5] 

2. Data 

We observed HD 167971 with the NRAO Very Large ArrajQ 
(VLA) on 2002 March 24 (at 3.6, 6 and 20 cm) and 2002 
September 11 (at 3.6, 6, 18 and 20 cm). All observations were 
centred near HD 167971, except the 2002 September 11 obser- 
vations at 18 and 20 cm, which were centred near HD 168112. 
The field covered at these wavelengths is large enough to con- 
tain HD 167971 as well. All observations consist of a single run 
on the target, preceded and followed by a phase calibrator. The 
details of the reduction are given in AppendixlAl 

In Paper I we described the Australia Telescope Compact 
ArrajQ (ATCA) observations which were collected on 2001 
October 11, centred on HD 168112. The 13 and 20 cm obser- 
vations have a field that is large enough to contain HD 167971. 
Table lA.ll gives more information about the reduction of this 
dataset and the fluxes we determined. 

We supplemented our own observations with archive data. 
The VLA archive contains a number of observations that are 
centred on, or close to, HD 167971. (The ATCA archive was 
also searched, but no additional data were found.) Many of these 
observations have not been published previously. To avoid intro- 
ducing systematic effects in the data reduction, we re-reduced all 

^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 

^ The Australia Telescope Compact Array is part of the Australia 
Telescope which is funded by the Commonwealth of Australia for op- 
eration as a National Facility managed by CSIRO. 



observations consistently. Details of the reduction and the mea- 
sured fluxes are given in Appendix lAl The 2, 3.6, 6 and 20 cm 
fluxes are plotted in Fig.[T] 

Fig. [1] confirms the non-thermal nature of the radio emis- 
sion from HD 167971. The fluxes increase with larger wave- 
length (with some exceptions at 20 cm) and the spectral index is 
therefore negative. The fluxes are also clearly variable and very 
high compared to what would be expected from thermal free- 
free emission. The stellar wind of the 08 supergiant dominates 
the free-free emission, with a thermal flux of only 0.1 mJy at 
6 cm based on typical wind parameters listed in Table [1^. 

3. Variability 

3. 1 . Possible instrumental effects 

A number of problems arose during the data reduction (see 
Appendix [All. One of those was that many observations have a 
phase calibrator at a rather large angular distance from the tar- 
get (9-11°), which forced us to apply selfcalibration to almost all 
observations. It is therefore important to check that the variabil- 
ity seen in Fig.[T]is not instrumental. First of all, the HD 167971 
flux should not be correlated with the flux derived for the phase 
calibrator The phase calibrator fluxes were determined from the 
AIPS task GETJY (except for the ATCA data, where we measured 
them on an image) and are listed in Table lA.ll We checked 
those calibrators that were used in three or more observations 
and found that the phase calibrator and the HD 167971 flux are 
not significantly correlated. 

Secondly, any flux changes of other targets on the same 
image should not be correlated with the flux changes of 
HD 167971. This test could only be applied at 6 cm (one other 
target) and 20 cm (three other targets). Again, when we tested 
this, no significant correlation was found. Furthermore, the range 
in variability in HD 167971 is substantially larger than that seen 
in the other targets. 

The VLA instrument can be used in a number of configu- 
rations that correspond to different spatial resolutions. We there- 
fore also need to check that the higher resolution observations do 
not partly resolve the target and therefore result in a lower flux. 
From VLBI observations, Phillips & Titus (1990) found that the 
linear size of the non-thermal region is at most 16 milh-arcsec 
(mas) at the time of their observation, which is well within the 
beamsize of the VLA (~ (X'S at 6 cm, in the highest resolution 
configuration). Furthermore, at 6 cm, the highest-resolution A 
configuration gives both the highest flux value and a very low 
one (see Table lA.ll ). At 20 cm, the A configuration gives the 
highest value, while the low values were obtained rather with 
the lower-resolution C configuration. The observed variability 
can therefore not be ascribed to our partly resolving HD 1 6797 1 . 

3.2. Long-term variability 

It is quite obvious from Fig. [T]that the fluxes show long-term 
variability. This is most clearly seen at 6 and 20 cm. The 6 cm 
fluxes vary between ~8 and ~ 18 mJy (or even higher, if we 
include the less reliable AC216 and AD219 observations) and 
the 20 cm fluxes vary between ~ 4 and ~ 25 mJy (over a time- 
scale of 5-10 years). 

The 6 cm dataset is the most complete, in the sense that it 
covers the ~ 20-yr time-range reasonably well. Judging by eye, 
a sine function seems quite appropriate to fit the long-term trend 
in these data (although this is based in part on the 1988 AC216 
and AD219 observations, which have large error bars). Fitting 
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Fig. 1. The 2, 3.6, 6 and 20 cm radio fluxes of HD 167971, plotted as a function of time. Not all 20 cm upper limits are included in 
this figure, only those which are lower than 40 mJy. Fluxes at 0.7, 1.3, 13 and 90 cm are not plotted, because only one or two data 
points are available. Note that there are differences of a few percent in wavelength for observations we classed in the same band 
(these differences are up to 20 % in the 20 cm band). Observations are labelled by their programme name. 



such a sine function, we found a ~ 20-yr period, with a maximum 
around 1986.2. Of course, a period that is so similar to the time- 
range covered by the observations is not really convincing: much 
longer periods could fit the data equally well. In that case one 
would have to assume that the fluxes are higher outside the time- 
range covered. 



The 20 cm data cannot be fitted with a 20-yr period sine func- 
tion. Either a much longer period would be needed, or the data 
do not follow a sine curve. From other colliding-wind binaries, 
such as WR 140 (Dougherty et al. |2005l l, we know that the flux 
variation at larger wavelengths does indeed show a sharper peak 
than at smaller wavelengths (see Sect. 14. 11 1. The behaviour of the 
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Fig. 2. Various subsets of radio observations are plotted as a 
function of orbital phase, with a period of f = 3.3213 days. 



HD 167971 20 cm fluxes seems consistent with this and there- 
fore does not contradict a 20-yr period. 

The 2 and 3.6 cm fluxes show some variability, but the range 
is small compared to the 6 and 20 cm range. The time coverage 
at 2 and 3.6 cm is not as good as for 6 cm, but is still good enough 
to suggest that this smaller variability is real and is not due to the 
limited sampling. 

The spectral index between 2 (or 3.6) and 6 cm is always 
negative, and remains between -0.3 and -0.6 most of the time. 
The 6-20 cm index is quite variable: it is -0.3 around 1985, 
is positive around 1990 (if we rely on the sine interpolation 
of the 6 cm data), is -0.0 in 1993-1994 and is -0.3 again 
in ~ 2002-2003. At certain times during the cycle (e.g. 1990- 
1994), the fluxes therefore no longer follow a power law, but 
show a turnover (from negative to positive spectral index) be- 
tween 6 and 20 cm. 



3.3. Short-term variability 

Superimposed on the long-term variability, there is also shorter- 
term variability with a smaller amplitude. Clear examples of this 
are seen in the 6 cm 1984-1985 observations (projects AA28, 
AA29 and AC116) and the 3.6 and 6 cm 1993-1994 observa- 
tions (projects AB671 and AL320). We will consider these two 
intervals separately: this has the advantage that we do not need 
to de-trend our data for the long-term variations. An important 
question that will be addressed is whether we can detect the 3.3- 
day period of the eclipsing binary in these data. 

We first of all note that these short-term variations are signifi- 
cant. This is most easily seen by determining the;^^^ from fitting a 
constant to the observations. We find a reduced;^^^ of 10.8 (1984- 
1985), 19.6 (1993-1994, 3.6 cm) and 14.0 (1993-1994, 6 cm). 



The error bars used in this calculation are from Table I A. ll They 
cover not only the root-mean-square (RMS) noise in the map, 
but also include the uncertainty in the absolute flux calibration 
and an estimate of some of the systematic errors (see Paper I). 
For observations with the smallest error bars, it is the absolute 
flux calibration that dominates. 

Fig.|2]shows that the 1984-1985 6 cm data show good orbital 
phase locking, suggesting that they indeed follow the 3.3-day 
period. However, with only 5 observations this good fit may not 
be very significant. A series of 100 Monte Carlo simulations in 
which the same fluxes were attributed random phases showed 
that about 30 % of those showed good phase locking (as judged 
by eye). 

For the 1993-1994 observations, neither the 3.6 nor 6 cm 
fluxes show good phase locking. Especially around phase 0.2- 
0.4 there is a large variation in flux, which is not compatible with 
the expected simple behaviour of the fluxes with orbital phase. If 
one would presume some of the observations around phase 0.2- 
0.4 to be in error and eliminate them, a further problem would 
arise: the phase of minimum flux would then shift considerably 
(by ~ 0.2) between the 1984-1985 and the (remaining) 1993- 
1994 observations. To verify if such a phase shift is possible, 
we checked the optical data from the Long-term Photometry of 
Variables campaign (LTPV - Manfroid et al. 119911 Sterken et 
al.ll993) and from Hipparcos (ESA |1997b . These data cover the 
time-range 1985-1993 and show a period that is slightly differ- 
ent from the Leitherer et al. ( 119871 ) one: we find 3.32157 days 
instead of 3.3213. This new value is closer to more recent de- 
terminations of the photometric period (3.321609 d, Mayer et 
al. igg?; 3.321634 d. Van Leeuwen & Van Genderen ,19971 l. 
Using even the most different periods (3.321634 d vs. 3.3213 d) 
results in a phase shift of only 0.1 over 9 years. This is too small 
to explain the required phase shift between the 1984-1985 and 
1993-1994 data. 

In summary, we conclude that the 3. 3 -day period of the 
eclipsing binary is not detected in the radio data. 

We also searched for other periods in the data. We de- 
trended the 3.6 and 6 cm observations by subtracting the best- 
fit sine curve with a 20-yr period. We then systematically tried 
periods between 3.3 d and 20 yr and evaluated the good- 
ness of phase locking with the minimum string-length method 
(Dworetskv l 19831 1. The most promising solutions were evaluated 
by eye, but no period was found that gave good phase locking for 
both the 3.6 and 6 cm data (in all cases, at least one data point 
was significantly discrepant). 



4. Interpretation 

4.1. Colliding-wind region inside the binary 

Here, we consider the possible colliding-wind region inside the 
binary. Because both binary components are of equal spectral 
type (Leitherer et al. I1987I I. the collision will occur at the mid- 
plane between the two stars, i.e. at 21 Rq for an 05 V -H 05 V 
binary (see Table[T^). If this colliding-wind region has an impor- 
tant contribution to the total flux, the 3. 3 -day binary period may 
be detectable in the fluxes. 

We did not find the 3.3-day period in the observations, how- 
ever. Two possible explanations can be suggested for this: either 
there is little or no intrinsic synchrotron emission, or there is too 
much free-free absorption for the synchrotron emission to es- 
cape. A lack of intrinsic emission could have various causes: the 
colliding-wind region could be too weak to provide enough en- 
ergy for the synchrotron emission, the Fermi acceleration mech- 
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Table 1. The top panel (a) shows typical parameters for an 05 V, 
an 08 V and an 08 I star, based on stars of similar spectral type 
found in Puis et al. ( I1996K Repolust et al. (2004) and Markova et 
al. (2004). We list the radius (Rt), mass (M*), terminal velocity 
(Voo) and mass-loss rate (M). The bottom panel (b) shows the or- 
bital parameters, both for the inner binary and for the third com- 
ponent assumed to be orbiting the inner binary. In both cases, 
we list two possibilities for the spectral types of the inner binary 
components. We list the semi-major axis (01,2) and the ampli- 
tude of the orbital motion {Ki^2), which were calculated from the 
mass listed in the top panel and an assumed inclination angle 
of 90°. The masses of an 08 V -H 08 V binary are rather low; 
they are not compatible with the order of 100 Mq for the total 
mass of the binary derived by Leitherer et al. ( 119871 ) from the 
observed velocity amplitude of ~ 300 km s ' . For the inner bi- 
nary the known period was used, while for the inner binary -1- 
08 I system, a 20-yr period was assumed. 
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Fig. 3. Contours of optical depth at 6 cm in the stellar wind of 
an 05 V star (left panel) and an 08 V star (right panel). The 
observer is situated to the left. The contour lines of r = 1, 10 
and 100 are shown. The two dots indicate the position of both 
components of the binary. 
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anism could be inefficient, or the inverse-Compton cooling could 
be too strong. We explore each of these possibilities in detail be- 
low. 

Even for such a close binary, the velocities of the material be- 
tween both components are supersonic and a coUiding-wind re- 
gion is formed (as X-ray observations of similar systems, such as 
the 3.4-day period HD 159176 confirm; De Becker et al. 2004). 
Radiative inhibition, however, considerably reduces the outflow 
velocities. Stevens & Pollock ( 119941 ) calculated the effect of ra- 
diative inhibition for a somewhat similar binary (HD 160652, 
06.5 V + 06.5 V, P - 6.14 d). Extrapolating the results of their 
Fig. 4 to our even closer system, we estimate that the outflow 
velocity at the shock is v ^ 500 km s"' . 

We first check that the collision between the winds is ener- 
getic enough to provide the synchrotron emission. In this order- 
of-magnitude analysis, we assume HD 167971 to be an 05 V 
H- 05 V binary with star and wind parameters based on litera- 
ture values for similar stars (see Table [TJ. The total kinetic en- 
ergy of both winds is Mv^ = 2.4 x lO"*^ erg s"'. Assuming the 
size of the emission region is ~7rroB (Eichler & Usov 1993), 
where roB - 21/?0, the fractional solid angle of the stellar 
winds impinging on the non-thermal emission region is ~ 0.23, 
and the kinetic luminosity available to the wind-collision region 
is 5.5 X lO-'^ erg s"'. We estimate a total radio luminosity of 
~ 1.6 X 10^" erg s"' at a distance of 2 kpc, which can easily be 
provided from the kinetic luminosity of the winds. 



For the Fermi acceleration to be efficient, the Alfven Mach 
number 



Ma = 



47rp 



B 



(1) 



should not be too small ( Wentzel ^19741 Jones & EUison IT99T] 
Sect. 4.1.1). B is the magnetic field at the shock and v and p are 
the pre-shock velocity and density. Assuming a value of Ma - 4 
and a pre-shock density of 7.0 x lO"'** g cm"^, we find that B 
should be 12 G. Assuming a Weber & Davis ( |1967| l type mag- 
netic field structure, this leads to a surface field that is small 
compared to the usual assumption of a 100 G surface field (com- 
patible with the detection limit, see Mathys 1999). 

We next explore the possibility that inverse-Compton (IC) 
cooling is too strong because of the large supply of UV photons 
available in this close binary. These photons IC scatter off the rel- 
ativistic electrons and thereby remove energy. The time needed 
for IC cooling to reduce the momentum of an electron from a 
high value pc down to p can be found by integrating Eq. (6) of 
Van Loo et al. (120051) : 



pcrjLt 



At : 



(2) 



where we assumed that p «: As we typically consider ra- 
dio fluxes in the 2-20 cm wavelength range, we determine the 
time required to cool down to a momentum that has maximum 
synchrotron emission just beyond 20 cm. Using the relation be- 
tween momentum and the frequency at which the maximum flux 
is emitted (Van Loo et al., their Eq. (19)), we find At ^ 600 s 
(for an assumed L, = IO^Lq). During that time, the electron has 
moved over a distance vAt » 0.5 Rq. 

Due to the IC cooling, the size of the synchrotron emitting 
region is therefore rather small. This makes it doubtful whether 
the colliding-wind region inside the binary can generate a signif- 
icant amount of synchrotron emission. Redoing the analysis for 
later-type components in the binary, or (super)giants, results in 
a similar size for the synchrotron emitting region (within a fac- 
tor of three). If the components are 05 supergiants, then we note 
that the Alfven number requirement results in a ~ 40 G magnetic 
field at the shock, which translates into a surface field that is in 
better agreement with an assumed surface field of 100 G. 

We finally explore the possibility that the free-free absorp- 
tion is so large that all of the synchrotron emission is absorbed 
at all phases in the orbit. In Fig. [3] we show the contours of op- 
tical depth (t) at 6 cm. To simplify the figure and discussion. 
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we consider the free-free effects of only a single star. The op- 
tical depth was calculated from the formulae given in Wright 
& Barlow ( |1975l l. Fig. [3] also shows the position of both com- 
ponents. For the purposes of this figure we have assumed that 
the inclination is exactly 90° (we know it is close to 90°, be- 
cause it is an eclipsing binary). The size of the colliding-wind 
region is somewhat larger than the binary separation and all of 
the synchrotron radiation will be found close to this region (see 
above). We see on Fig.|3]that the synchrotron-emitting region is 
well within the t = 10 contour surface and would therefore be 
undetectable (If the components are both giants or supergiants 
(Davidge & Forbes 1988) rather than main-sequence stars, then 
the opacity effects discussed here will be even stronger) 

Another possibility is the intrinsic synchrotron emission 
from the wind-collision region in the binary is absorbed by free- 
free opacity of the 08 supergiant wind. However, as the relative 
positions of the supergiant and the binary system are unknown 
this remains an open question. 

4.2. Colliding-wind region between the supergiant and the 
binary 

The major part of the observed non-thermal emission of 
HD 167971 is not due to the colliding-wind region inside the 
binary, because (i) we do not detect the 3.3-day period, (ii) the 
intrinsic synchrotron flux is estimated to be small, and (iii) the 
free-free absorption is large. We can therefore conclude that the 
main source of the observed non-thermal radio emission must be 
the colliding-wind region between the 08 supergiant and the bi- 
nary. We recall that from spectroscopic observations (Leitherer 
et al. 1 19871 1 it is not clear if the supergiant and the binary are 
gravitationally bound. The supergiant could be coincidentally in 
the same line of sight and might even be at a (somewhat) dif- 
ferent distance. However, we will assume in our analysis that 
flie supergiant and the eclipsing binary do form a gravitationally 
bound system. This allows us to attribute the long-term variabil- 
ity to the relative motion of the binary and the supergiant. 

The orbital period is then ~20 years, or possibly longer 
(Sect. 13.2b . With the typical stellar parameters for an 08 super- 
giant and the eclipsing binary, a 20-year period corresponds to a 
separation 01+02 > 6700 Rq (see Table [TJ)), provided we as- 
sume a nearly-circular orbit. This is well beyond the ~ 1600 Rq 
radius at which t = 1 for 6 cm in the supergiant wind. Most of 
the time we therefore see the intrinsic synchrotron emission at 
this wavelength. 

At 20 cm, the r = 1 radius is much larger (~ 3500 Rq). A 
much longer part of the orbit will therefore be spent inside the 
20 cm radio photosphere (except for inclinations near 0°), ex- 
plaining why the maximum at 20 cm is much more peaked. At 
wavelengths shorter than 6 cm, the radio photosphere is much 
smaller and the colliding-wind region is therefore outside the 
radio photosphere nearly all the time, explaining the lack of sub- 
stantial variation at 2 and 3.6 cm. The observed changes in the 
6-20 cm spectral index and the turnover between 6 and 20 cm 
during part of the period (Sect. I3.2l i can thus be ascribed to 
these opacity effects. An alternative explanation for the turnover 



One should be careful with this conclusion. In the WR 140 sys- 
tem, for instance, the colliding-wind region is well within the free-free 
radius of the Wolf-Rayet star during the whole orbit. The O star wind 
blows a low-density region in the Wolf-Rayet wind and it is through this 
region that we can see the synchrotron emission, at least part of the time 
(Williams et al. 1 1990b . The winds of the HD 167971 binary are nearly 
equal, and therefore no such complications arise here. 



would be the Razin effect. Due to this effect, the synchrotron 
emission deviates from a power law at longer wavelengths and 
the long-wavelength spectral index can become zero or even 
positive. The strength of the Razin effect depends on the (lo- 
cal) magnetic field and electron density. In an eccentric orbit, 
both quantities would vary with orbital phase, and therefore, so 
would the spectral index. More sophisticated models would be 
required to derive quantitative orbital information from this al- 
ternative explanation of the turnover 

Eichler & Usov ( 119931) made a simple model for the syn- 
chrotron radiation emitted in a colliding-wind binary. When we 
apply their analysis to the present data, we find an intrinsic syn- 
chrotron flux (i.e. without free-free absorption) of ~ 25 mJy at 
6 cm, which is in surprisingly good agreement with the observed 
values. The Eichler & Usov model does contain a number of un- 
constrained parameters for which we adopted their default val- 
ues. Hence, the agreement should only be interpreted as showing 
the plausibility of our model. 

We next try to explain the variability of the flux and the 
turnover between 6 and 20 cm (which occurs during part of 
the orbit) by ascribing them to the changing free-free absorp- 
tion due to the orbital motion of the synchrotron-emitting source 
and the 08 supergiant. We used a simple model, where we as- 
sumed the intrinsic synchrotron emission to come from a point 
source in a 20-yr circular orbit around the supergiant (with stel- 
lar parameters as in Table [TK). Simulations show that we can 
reproduce the 6 cm radio lightcurve quite well, provided we use 
a low-inclination orbit (otherwise the flux would be completely 
absorbed). However, free-free absorption increases considerably 
from 6 to 20 cm, and we cannot match the 6 and 20 cm data with 
the same model. This is most probably due to our assumption 
that the synchrotron-emitting region is a point source. In real- 
ity, the synchrotron-emitting region will be extended and deriv- 
ing quantitative information about the opacity from the flux or 
the turnover will therefore require more sophisticated models, 
which are beyond the scope of this paper 

Further information on the size of the synchrotron-emitting 
region comes from the Phillips & Titus ( 1990) observations. 
They found a size of 16 mas, which corresponds to a diame- 
ter of ~6500 Rq at 18 cm. This is comparable to the separa- 
tion between the binary and the O supergiant. The size of the 
colliding-wind region makes it difficult to put constraints on the 
orbital inclination. If the inclination of the orbit is close to 90°, 
a point-like source of synchrotron emission would be eclipsed at 
6 cm during ~ 2 years of the 20-year orbit, as it passes behind 
the free-free radius of the supergiant; but this is contradicted by 
the observations. However, for a larger-sized colliding-wind re- 
gion, this argument cannot be used, as there is probably always 
some part of it that is not eclipsed. No constraints on the orbital 
inclination can therefore be derived. 

The shorter-term variability could be attributed to stochas- 
tic variations due to clumping in the stellar winds, as suggested 
for WR 146 and WR 147 (Setia Gunawan et al. 2000, 2001). 
These affect the synchrotron emitting region itself (as the den- 
sity and/or velocity of the material entering the colliding-wind 
region changes) and affect the free-free absorption. The AB67 1 
series represents a typical example of such flux variations. The 
6 cm fluxes range between 10 mJy and 12 mJy. We can derive 
an estimate of how much the density needs to change to explain 
this flux variation through additional free-free absorption. If, for 
simplicity, we assume the synchrotron emitting region to be a 
point source, we can use the Wright & Barlow ( |I975| l formu- 
lae to determine how much optical depth has to be added to 
change the 12 mJy flux into 10 mJy. When the free-free opti- 
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cal depth along the Une-of-sight is 1, a 10 % increase in density 
suffices. For larger-sized sources, the density contrasts will need 
to be considerably larger to offset the cancellations that occur in 
a stochastically clumped wind. 

The model proposed here is also compatible with the X-ray 
data, studied by De Becker et al. (I2005bl) . From XMM-Newton 
observations, they found an excess of a factor of 4, compared to 
the canonical Lx vs. Lboi relation for O stars. The emission is 
probably thermal and shows the presence of a high-temperature 
plasma component (~ 20 MK). The colliding-wind region within 
the eclipsing binary is unlikely to produce such high plasma 
temperatures as the winds of the two 05-08 V stars collide be- 
fore they reach their terminal velocities, thereby producing softer 
thermal X-ray emission. The high plasma temperature should 
therefore be ascribed to material heated in the colliding-wind re- 
gion between the eclipsing binary and the 08 supergiant. The 
existence of this interaction shows that the supergiant must be 
at approximately the same distance as the binary, and is not an 
accidental line-of-sight object. 

The two ROS AT data points plotted on Fig. 5 of De Becker et 
al. ('2005b) show X-ray flux decreasing between 1993 and 1995. 
This decline closely follows the radio emission during the same 
period. We attribute this to the intrinsic X-ray and synchrotron 
emission varying in a correlated way, as they both depend on the 
shock strength and the absorption in the stellar wind material. 
However, De Becker et al. caution that the observed variability 
might be attributed to systematic effects between XMM-Newton 
an d ROSA T data, as seen in Cyg OB2 No. 8A (De Becker et 
al. l2005al l. 

5. Conclusions 

In this paper we analysed the available radio data on HD 16797 1 
from the VLA and ATCA. The negative spectral index confirms 
the non-thermal nature of the radio emission. The fluxes are also 
very high (compared to thermal) and they are variable. 

At the core of the HD 167971 triple system is a 3.3-day pe- 
riod eclipsing binary. No modulation of the radio fluxes with 
this period was detected. We can explain this by a rather low 
intrinsic synchrotron flux (due to a low acceleration efficiency 
and/or strong inverse-Compton cooling), or by the large amount 
of free-free absorption in the stellar wind material. We there- 
fore conclude that most of the non-thermal emission must come 
from the colliding-wind region between the OS supergiant and 
the combined winds of the binary. 

The hypothesis that the 08 supergiant and the binary have 
a colliding-wind interaction provides a plausible explanation for 
many of the characteristics of the radio and X-ray observations. 
If, furthermore, the system is gravitationally bound, the observed 
variability can be ascribed to the orbital motion. The proposed 
orbital period of ~ 20 years is large enough to be well beyond the 
T = 1 free-free radius. Some absorption still occurs, explaining 
the long-term variability of the radio fluxes. The high X-ray tem- 
perature plasma component found by De Becker et al. ( 2005b) 
can also be ascribed to this colliding-wind region. The problem 
of fitting the 6 and 20 cm radio fluxes simultaneously with a 
point-source model, and the Phillips & Titus ( 1990) resolved ob- 
servations suggest a colliding-wind region of substantial size. 

While this model provides a good explanation for the ob- 
servations, independent verification still needs to be sought. 
Verification of our HD 167971 model can most easily be ob- 
tained from high-resolution, high S/N optical spectra. These 
should bring definite evidence that the triple system is indeed 
gravitationally bound and should allow us to obtain better stellar 



and orbital parameters for the binary. Once these data are avail- 
able, quantitative models for the X-ray and non-thermal radio 
emission can be developed. 
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Appendix A: Data reduction 

The data reduction was done using the Astronomical Image 
Processing System (AIPS), developed by the NRAO. The dif- 
ferent reduction steps consist of assigning the standard fluxes 
to the flux calibrators, calibration, imaging and cleaning. These 
steps have been described in detail in the Appendix to Paper I. 
On some observations made with low spatial resolution. Galactic 
background structure can be seen. This background can easily be 
filtered out by using the properties of the Fourier transform that 
is an essential part of radio interferometry: we simply drop those 
visibilities that were measured on the shortest baselines and 
thereby eliminate the large-scale background structure. In doing 
so, we took care not to remove more data than strictly needed. 
Technical information on the reduction is listed in Table I A. f| 

In a considerable number of cases, the distance of the phase 
calibrator to HD 167971 is large (9-11°, see Table lAlTl. The an- 
tenna gains could therefore be significantly different between the 
phase calibrator and target scans. We therefore systematically 
applied selfcalibration to those observations that have detections. 
We did one round of phase-only selfcalibration (a second round 
does not make any significant difference). For the time-interval 
over which to integrate, we took the shortest time-interval that 
gave smooth phase solutions (phase scatter of less ffian ~ 20°): 
this is usually 3 or 5 min. Observations with phase scatter sig- 
nificantly larger than 20° have been listed in the "notes" column 
to Table lA.ll No selfcalibration was applied to the AC308 ob- 
servations (because of their short on-target time) and the ATCA 
data. 

On the selfcalibrated images, the fluxes were measured by 
fitting an elliptical Gaussian to the source. We checked that the 
resulting fluxes are in good agreement with the maximum inten- 
sity (this is consistent with HD 167971 being a point source on 
all our images). The fluxes are listed in Table lA.ll The necessity 
for using selfcalibration is clearly shown by a number of obser- 
vations that have substantially larger fluxes for HD 167971 on 
the selfcalibrated images than on the non-selfcalibrated images. 
The largest effect is seen in the AB671 (1993-02-01) 3.6 cm ob- 
servation, where the selfcalibrated flux is 6.2 + 1 mJy, while the 
non-selfcalibrated one is only 2.6 mJy. 

The error bars listed in this table cover not only ffie root- 
mean-square (RMS) noise in the map, but also include the un- 
certainty in the absolute flux calibration (2 % for the 3.6, 6, 20 
and 90 cm observations and 5 % for the other wavelengths) and 
an estimate of some of the systematic errors (see Paper I). For 
those observations with relatively small error bars, it is the un- 
certainty in the absolute flux calibration that dominates. In those 
cases where HD 167971 is off-centre, its flux was corrected for 
primary beam attenuation and beamwidth smearing. 

For the observations where HD 167971 was not detected, 
we assign an upper limit of 3 times the RMS, where the RMS 
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is measured in a small box around the source. If the source is 
off-centre, this upper limit is corrected for the primary beam at- 
tenuation and beamwidth smearing. A number of 20 cm observa- 
tions containing HD 167971 are of such low quality that they do 
not even provide a significant upper limit: BAUD (1980-07-12 
and 1982-02-28), FIX (1982-05-22) and AT143 (1993-06-11). 
A number of 90 cm observations were also rejected for the same 
reason. 

In some 20 cm observations HD 167971 is within the pri- 
mary beam on two images centred on different targets. In that 
case we only list the one closest to HD 167971 because this 
gives the smallest error bar. In most cases the flux of the other 
image is compatible with the flux listed. The exceptions are: 
AC116 (1984-DEC-21), where the off'-centre measurement of 
HD 167971 is a factor of two lower than the on-centre mea- 
surement and BPl, where the off'-centre value is a factor of 
two higher than the on-centre value. A smaller eff'ect is seen 
for VP51, where the oflF-centre value is 30 % higher than the 
on-centre value. Simulations show that the efifect in the high- 
resolution ACl 16 observation is due to noise in the map. For the 
low -resolution BPl and VP51, we checked that decreasing the 
background even further (by throwing away more short base- 
lines) leads to a much better agreement. Some caution must 
therefore be used in off'-centre fluxes, such as the low-resolution 
AC308 observations. (Note that the values Hsted for ACl 16, BPl 
and VP51 in Table lA.fl do not suffer from this problem because 
they are based on the on-centre measurements.) 

A number of observations have used two ffux calibrators. 
We checked that using either calibrator results in the same ffux 
(within the error bars). This gave good agreement in most cases, 
except AL320 (1994-FEB-18) 3.6 cm, where there is a 32 % 
difference. Checking other observations made around the same 
time shows that the (uncalibrated) ffuxes of 3C48 are anoma- 
lously low for the observation we are interested in, while the 
(uncalibrated) 3C286 ffuxes remain the same over all the obser- 
vations considered. The calibration of this observation is there- 
fore made on 3C286. A similar problem exists for AL320 (1994- 
FEB-18) 6 cm. Applying the same technique shows that, in this 
case, 3C48 is the more reliable calibrator. We also note that for 
the TSTOB observation we used the ffux calibrator from another 
programme, nearby in time (see Paper I). 

The reduction of the ATCA observations C978 was detailed 
in Paper I. These observations are centred on HD 168 1 12, but the 
1 3 and 20 cm primary beams are large enough that they also con- 
tain HD 167971. Because these data were collected in a number 
of spectral channels, we can apply multi-frequency synthesis: 
this results in images that are not beamwidth smeared. 

A comparison between our ffuxes and those that have been 
published in the literature shows that our error bars are usually 
larger, because we include the ffux calibration errors and some 
estimate of systematic effects. Our 2 cm determinations are sys- 
tematically higher than the literature values. This is because of 
the improvement due to the selfcalibration reduction technique. 
A similar, but smaller, eff'ect is seen at 6 cm. For the 20 cm 
observations of Bieging et al. d 19891 1 we obtain lower ffuxes. 
The VLA Calibrator Manual instructions (Perley & Tavlor |2003T l 
for this specific combination of wavelength, ffux calibrator and 
VLA configuration show that the results obtained should be re- 
duced by 6 %. We speculate that Bieging et al. did not apply this 
6 % reduction; doing so results in an acceptable agreement be- 
tween their ffuxes and ours. Our upper limit to the AR328 0.7 cm 
observation (4.5 mJy) is much higher than the Contreras et al. 
( 1 19961 1 value of 1.72 mJy, We note that this is a 19 min observa- 
tion using only 10 antennas, which gives a theoretical RMS of 



0.7 mJy/beam. The 3 sigma upper limit of 1.72 mJy claimed by 
Contreras et al. is therefore too optimistic, and we have greater 
confidence in our result. 
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Table A.l. Reduction of VLA and ATCA data. Programme C978 is an ATCA observation, all others are VLA observations. Column (1) gives the 
programme name, (2) the date of the observation, (3) the phase calibrator (J2000 coordinates), (4) the flux of the phase calibrator + the rms (in 
Jy), (5) the distance of HD 167971 to the phase calibrator (degrees), (6) the integration time (in minutes) on the source, (7) the number of antennas 
that gave a usable signal, (8) the configuration the VLA was in at the time of the observation, (9) the RMS in the centre of the image (not listed in 
case on an upper limit), (10) the measured flux (in mjy) and (11) refers to the notes. Many of the VLA observations were made in two sidebands, 
each of which has a bandwidth of 50 MHz; the exceptions are noted in column (1 1). Upper limits are 3 x the RMS. Numbers between brackets in 
the notes column give references for those observations that have already been published in the literature. 
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1 OOO 1 A/1 

iyyy-iz-u4 


1822-096 


1.37 ± 0.02 


2.8 


121 


19 


B 


A Al 


1 A A _i_ A 

lU.U + U.Z 


A 


T>T> 1 1 A 
DD 1 10 


onnn aa oa 
ZUUU-UO-ZO 


1822-096 


1.340 ± 0.003 


2.8 


117 


26 


DnC 


A Al 


O -u A 

y.z ± U.Z 




Ad lU4o 


onno ni ia 

ZUUZ-Uj-Z4 


1832-105 


1.67 ±0.04 


3.9 


8 


26 


A 


A AO 

U.U / 


1 1 _u A /I 
IZ.J ± u.4 




Ad 1 UO J 


onno no 1 1 
zuuz-uy- 1 1 


1832-105 


1.29 ±0.02 


3.9 


7 


26 


CnB 


A 1 1 

U.l 1 


C Q -1- A 1 

o.y ± U.J 




6 cm 

A A 

AAZo 


lQC/1 ni AO 

iy?54-U:>-Uy 


1733-130 


5.74 ± 0.01 


11.0 


19 


26 


CnB 


A AQ 
U.Ub 


1 c Q , n /I 
Ij.o ± U.4 


/■ 1 ^ 
(1) 


A A OO 
AAZy 


1QQ/1 n/1 A/1 
iyo4-U4-U4 


1733-130 


5.69 ± 0.03 


11.0 


11 


27 


C 


A 1 1 

U.l 1 


1 ^ A -u A 1 

Ij.U ± U.J 


/■i ^ 
(1) 


Ar^l 1 A 
AL.1 10 


1 QQA 1 1 oo 

iyo4-i 1-Z / 


1733-130 


5.02 ± 0.02 


11.0 


21 


24 


A 


A AQ 

u.uy 


1 O 1 _u A /I 
1 / .I ± U.4 


(1) 


A /"'I 1 A 

AL-l 10 


1 OQ/1 1 O 01 

iyo4-iz-zi 


1733-130 


5.05 + 0.02 


11.0 


41 


26 


A 


A AA 

U.UO 


1 A < _i_ A 1 
10. J + U.J 


(1) 


Ad 10 


1 QQ< AO 1 A 

iy<5j-uz-io 


1733-130 


5.312 + 0.015 


11.0 


19 


24 


A 


A AQ 

U.uy 


1 /I O _L A 1 
14. / + U.J 


^'^ \ 
W 


AUZIO 


1 QQQ AO oo 

iyoc5-UZ-Z/ 


1811-209 


0.303 ± 0.002 


8.8 


15 


25 


CnB 


A AO 

U.U/ 


OA _i_ 1 A 

ZO + lU 


D,rrl 


A T\1 1 O 

Auziy 


1 QOQ n/1 1 ^ 

iyob-U4-i J 


1804+010 


1.149 ±0.004 


13.7 


17 


26 


CD 


A A/1 
U.U4 


1/1 -u A 
j4 ± 


R DT4 


A RAO 1 
Add / 1 


1 QQQ ni 0/1 

iyyj-ui-z4 


1811-209 


0.324 ± 0.002 


8.8 


19 


27 


A 


A A^ 
U.U J 


1 A O _u A 
lU. / ± U.Z 




A T3 AO 1 
Add / 1 


1 QQQ ni OQ 

lyyj-ui-zy 


1811-209 


0.318 ±0.001 


8.8 


19 


26 


BnA 


A AO 
U.U / 


11 O -1- A 1 
1 1. / ± U.J 




A RAO 1 

Add / 1 


1 QQQ no Al 
lyyj-uz-ui 


1811-209 


0.314 ±0.001 


8.8 


13 


27 


BnA 


A AO 
U.U / 


1 A /I -1- A 
1U.4 ± U.Z 




ARA7 1 
ADD / 1 


1 QQ1 AO AA 

lyyj-uz-uo 


1811-209 


0.322 ± 0.001 


8.8 


19 


27 


BnA 


A AA 
U.UO 


Q O _i_ A 
y. / ± U.Z 




A RATI 
AdO / 1 


1 Q01 AO 1 A 

iyyj-uz-14 


1811-209 


0.319 ±0.001 


8.8 


10 


27 


BnA 


A AC 

U.Uo 


10 1 _u A 1 
IZ.l ± U.J 




AT ion 


1 QQA AO 1 8 

iyy4-uz-io 


1733-130 


5.09 + 0.06 


11.0 


9 


26 


AC 


A AQ 

U.uy 


Q O _u A 1 
y. / ± U.J 


rxl 


AT ion 


1 QQA A< 1 A 

iyy4-uj-iu 


1733-130 


4.43 + 0.02 


11.0 


10 


26 


AB 


A AQ 
U.UtS 


1 A 1 _i. A 
lU.J ± U.Z 




AKjZo 


1 QO^ A4 oo 

iyyj-u4-z / 


1733-130 


5.04 ± 0.02 


11.0 


11 


16 


D 


A OA 
U.ZU 


S A -u A 4 
o.O ± U.4 


w) 


AW J 1 J 


I QQQ AA AC 

lyyy-uo-uo 


I OQ T 1 AC 

1832-105 


1.200 ±0.015 


3.9 


10 


20 


AD 


A 1 A 
U. 10 


1 A ^ -u f\ A 
lU.j ± U.4 




RR 1 1 A 
DD 1 10 


OAAA AA OA 
ZUUU-UO-ZO 


1822-096 


2.25 ± 0.03 


2.8 


93 


26 


DnC 


A A^ 

U.Uj 


11 ^ _i_ A 1 
1 1 . J ± U.J 


V 
A 


AB1048 


2002-03-24 


1832-105 


1.32 ±0.01 


3.9 


6 


26 


A 


0.09 


15.1 ± 0.3 




AB1065 


2002-09-11 


1832-105 


1.26 ± 0.07 


3.9 


5 


26 


CnB 


0.14 


13.4 + 0.3 




13 cm 






















C978 


2001-10-11 


1832-105 


1.042 ± 0.002 


3.9 120 


6 


EW352 


0.26 


11+2 


B,X 


20 cm 






















AC116 


1984-11-27 


1733-130 


5.55 ± 0.03 


11.0 


28 


25 


A 


0.12 


24.3 ± 0.6 


PH,(1) 


AC116 


1984-12-21 


1733-130 


5.90 ± 0.04 


11.0 


10 


27 


A 


0.17 


24.9 + 0.6 


(1) 


AG163 


1984-12-24 


1743-038 


1.281 ±0.015 


11.9 


1 


27 


A 




< 100 


B,2x25 


AC116 


1985-02-16 


1733-130 


6.36 ± 0.03 


11.0 


15 


25 


A 


0.12 


22.0 + 0.5 


(1) 


AB531 


1989-05-01 


1834-126 


0.518 ±0.001 


3.9 


2 


26 


B 




< 100 


B,lx50 


VP51 


1989-06-19 


1733-130 


5.200 


11.0 


59 


25 


C 


0.24 


6.9 + 0.3 


PH,X,C 


AG290 


1989-07-20 


1833-210 


13.02 ± 0.03 


9.6 





27 


CD 




< 80 


B 


VP91 


1989-11-10 


1733-130 


6.2 ± 0.2 


11.0 


56 


27 


D 


1.5 


4.2 + 2 


PH.X.(2) 


BPl 


1992-05-30 


1911-201 


3.0 ± 0.2 


15.0 


91 


25 


CD 


0.8 


5.6 ± 3 


PH 



10 

Table A.l. continued. 



R. Blomme et al.: Non-thermal radio emission from 0-type stars. II 



(1) 


(2) 


(3) 




(4) 




(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


nrosr 


date 




phase calibrator 




inter 


no. 


config 


RMS 


flux 


notes 






name 


flux 


dist. 


time 


ants. 




(mJv) 


(mJy) 




TSTOB 


1993-07-05 


1822-096 


4.3 


+ 0.2 


2.i 


13 


27 


c 


0.3 


12.9 ± 0.9 




AS534 


1994-05-21 




no A 


5.39 ± 


0.04 


2.8 


4 


21 


BnA 




< 70 


B, 1x3. 125 


AS534 


1994-05-26 


1822- 


096 


4.99 ± 


0.03 


2.8 


5 


25 


BnA 




< 34 


B,X,lx3.125 


AS534 


1994-06-01 


1822- 


096 


5.15 ± 


0.02 


2.8 


4 


25 


BnA 




< 150 


B,X,lx3.125 


AC308 


1996-06-08 


1833- 


210 


iU. / o 3: 


n CO 


9.6 


1 


27 


DnC 




< 190 


B 


AC308 


1996-06-09 


1833- 


210 


10.80 ± 


0.02 


9.6 


1 


27 


DnC 


2.0 


12 ± 4 


B 


AC308 


1996-06-16 


1833- 


210 


10.78 ± 


0.02 


9.6 




27 


DnC 




< 70 


B 


AC308 


1997-10-16 


1833- 


210 


10.07 ± 


0.02 


9.6 


1 


27 


DnC 


2.0 


12 ±6 


B 


C978 


2001-10-11 


1832- 


105 


0.929 ± 


0.001 


3.9 


120 


6 


EW352 


0.42 


14 ±2 


B,X 


AB1048 


2002-03-24 


1834- 


126 


0.534 ± 


0.001 


3.9 


3 


26 


A 


0.20 


20.9 ± 0.6 




AB1065 


2002-09-11 


1834- 


126 


0.468 ± 


0.001 


3.9 


5 


26 


CnB 


6 


29 ±23 


B,I,18cm 


AB1065 


2002-09-11 


1834- 


126 


0.529 ± 


0.001 


3.9 


5 


21 


CnB 


3 


27 + 7 


B,20cm 


90 cm 


























AH299 


1988-06-21 


1829+487 


42.6 ± 


0.5 


61.0 


47 


21 


DnC 




< 270 


B,2x3.125 


AH299 


1989-05-28 


1829+487 


41.7 ± 


0.5 


61.0 


57 


25 


CnB 




< 80 


B,2x3.125 



Notes: 
B 

C 

I 

PH 
X 



(1) 
(2) 
(3) 



Observation is not centred on HD 167971. The flux has been corrected for that. 
No correction for beamwidth smearing needs to be applied to C978. 

for VP51, no flux calibrator was available, so we calibrated the flux on the phase calibrator, 
to which we assigned the value listed in Perley & Taylor ( I2003t . 
Interference is high in the AB 1065 18 cm image. 

The phases in the selfcalibration gain solutions show scatter that is significantly larger than 20°. 
The phase calibrator is of low quality (AB671, C978), is not listed in the VLA Calibrator Manual 
(Perley & Tavlor ,2003i AW515 at 1.3 cm, BBl 16), or the recommended constraints on the 
uv-range could not be applied (AR328 at 3.6 cm). 
Bieging et al. OlM- 
Phillips & Titus (1 99 0]l. 
Contreras et al. (il996b. 
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